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Abstract 

The scaling of charged hadron fragmentation functions to the Tsallis distribution for the momentum fraction 0.01 
x J$ 0.2 is presented for various e + e~ collision energies. A possible microcanonical generalisation of the Tsallis 
distribution is proposed, which gives good agreement with measured data up to x « 1. The proposal is based on 
superstatistics and a Koba- Nielsen -Olesen (KNO) like scaling of multiplicity distributions in e + e~ experiments. 
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1. Introduction 

The question: 'why do thermal models work in high 
energy collisions' may be answered on a purely mathe- 
matical basis. Thermal models are based on the canoni- 
cal distribution e.g. the Boltzmann-Gibbs (BG) (or gen- 
eralisations of it |23l|20ll34ll ). The canonical BG distri- 
butions follows directly from the Large Deviations The- 
orem (LDT) [ 1 , 2 , 3 ] if the ensemble consists of a large 
number of independent and identically distributed par- 
ticles with fixed total energy. This theorem, however, 
does not deal with the reason behind the identical dis- 
tribution of these particles: whether it is due to the fact 
that they are thermal particles, or that they are an en- 
semble of identically produced particles. In any case, 
their total energy is conserved, either during their pro- 
duction, or during their thermalisation. This is probably 
the reason why Hagedorn's model on thermal hadron 
production [4] does not only work for heavy -ion col- 
lisions (for a few examples see Refs. [4|-[7|, and lfT2l ). 
but for e + e~ annihilations as well, which are conjectured 
to be non-equilibrium phenomena. 

However, the BG statistics describes uncorrelated and 
non-interacting particles, and neither of these features 
are characteristics of particles stemming from high en- 
ergy collisions. Apparently, the BG distribution fails 
to describe transverse momentum (pj) spectra for in- 
termediate and high momenta (p T gl 3 GeV/c), while 
the cut power-law Tsallis distribution (TS) gives a good 
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fit in a wide measured range of proton-proton (pp), 
proton-antiproton (pp), nucleus-nucleus (AA) [8|-[20|, 
and e + e~ ET\ l22l collisions. Furthermore, charged 
hadron fragmentation functions also scale to the TS dis- 
tribution for 0.01 55 x 0.2, as shown in Section [5] 
where x stands for the momentum fraction of the pro- 
duced hadron, x = 2p/ -\fs, with -\fs being the center of 
mass energy. 

Though numerous proposed explanations exist for the 
emergence of TS statistics in equilibrium and out of 
equilibrium systems (for reviews see Refs. [20, 23 1), the 
reason for its occurence in high energy physics has not 
yet been clarified. Moreover, the TS distribution, being 
a generalisation of the canonical BG, has a non-compact 
support, and thus cannot be valid close to x ~ 1, where 
it should have a cut due to the finiteness of the collision 
energy. Therefore, in order to achieve good agreement 
with data for large x, microcanonical effects, due to the 
finite available phase space, have to be taken into ac- 
count, as it has also been proposed in [24|-[29|. In this 
paper, a possible microcanonical generalisation of the 
TS distribution is proposed, which gives a good fit to 
data on fragmentation functions measured in e + <?~ colli- 
sions for 0.01 J$ x ;$ 1. 

This paper is organised as follows: in Section [2] it 
is shown that multiplicity distributions, similar to the 
Koba - Nielsen - Olesen (KNO) scaling ones reported in 
Ref. IT331 . can lead to a momentum distribution of the 
canonical TS type, or to a generalisation of it, if the 
momentum distribution is BG or microcanonical for a 
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fixed multiplicity. Section [3] contains fits of multiplic- 
ity distributions and fragmentation functions, used and 
obtained in Section[2] to e + <?~ data measured at various 
collision energies. Section|4]contains the concluding re- 
marks. 



2. Canonical and microcanonical Tsallis statistics 
from multiplicity fluctuations 

The fragmentation function measured in e + e~ experi- 
ments is the average of charged hadron yields produced 
in millions of events. For collision energies V* ^ M z 
(Mz is the mass of the Z boson), more than 90 % of 
the events are 2-jet events (a q and an q jet with the 
same energy, E = yfs/2). Hadrons in a q jet have a 
very narrow distribution arround the jet axis in momen- 
tum space, therefore they may be considered as a one- 
dimensional ensemble. Neglecting the mass of the pro- 
duced hadrons, energy and momentum conservation in 
one dimension has the same form: 2 e, = £ (note that 
all hadrons go in the same direction, and e,- = \p/\). If we 
prescribe nothing else than energy-momentum conser- 
vation during the hadronisation, the produced particles 
become a microcanonical ensemble in one dimension. 
In this section we show that for massless particles, 

a) if the momentum distribution in events with fixed 

multiplicity is BG, and the multiplicity has 
Gamma-distribution, the average momentum dis- 
tribution is TS-like; 

b) if the momentum distribution in events with fixed 

multiplicity is microcanonical, and the shifted mul- 
tiplicity (N - No) has Gamma-distribution, the ap- 
proximate average momentum distribution can be 
considered as a possible microcanonical generali- 
sation of the TS distribution. 

In case, a, the normalised one-particle distribution in 
events with multiplicity N, is 



Me) = A c e-^ 



(1) 



where A c - ff^j ' ( k D T(D) ) and k D — J dQ. p follow from 
the conditions 



E 
N 



= J" d£l p J" dpp D - l f N (e), 
= J d£l p f dpp D - 1 ef N {e). 



(2) 



ho is the angular part of the momentum space integral, 
D is the dimension of the phase space, T(x) is the Euler- 
Gamma function. From Eq. it follows that the in- 



(3) 



verse temperature in each event is proportional to the 
multiplicity: 

DN 
~E~' 

The multiplicity has Gamma distribution, which is the 
p = 1 case of Eq. (5j in Ref. [33]: 

p(N) =A m N a ~ x e~P N , (4) 

with fit parameters A m ,a,/3. In this case, the average 
momentum distribution is TS: 



kd,e 



a+D+l 



(5) 



The discrete sum has been approximated by a con- 
tinuous integral, x = e/E, the energy of the pro- 
duced hadron scaled by the maximal acquirable en- 
ergy in a 2-jet event E = yfs/2, and kd,e — 
A m D D /3 a+D+1 /(E D k D T(D)). Note that from Eqs. ^} 
and (|4]) it follows that the inverse temperature, fix, is 
also Gamma-distributed. This can also lead to TS mo- 
mentum distribution (see Refs. |[T8l[T9ll ). 

In case b, the normalised one-particle distribution in 
events with multiplicity N is 

Me)=A mc a-x) m - 1) -\ (6) 

where A mc = ( D fJ!.iy^D/(k D E D ) follows from the 
condition 



1 = J" dn p J dpp D - x f N {e). 



(7) 



Eq. ([6]) follows from the microcanonical momentum- 
space volume at fixed energy and multiplicity, 



N\ 



(8) 



from which the one-particle distribution is obtained as 

n N -i(E - e) 



Me) . 



(9) 



At this point, it is exploited that a shift in the multiplicity 
(N — > N - Nq) may be made without violating the KNO 
scaling (see Eq. (2) in 11331 ). and the averaging is done 
over the multiplicity distribution, 

p(N) = A m (N - Nof^e-^-^. (10) 

The resulting momentum distribution is a possible mi- 
crocanonical generalisation of the TS: 
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The discrete sum has again been replaced by a contin- 
uous integral, and the factorials in Eq. (|6]l have been 
approximated by the next-to-leading order Stirling- 
formula. The shift in the multiplicity has been chosen to 
be A^o = 1+2/D. This choice of No may seem rather ad 
hoc, but fittings shown in Section[3](Figs. T3fT6 1 justify 
it. 

A similar formula (Eq. 



1 1 1 without the 1 — x factor 



in the numerator) arises in the canonical approach for 
a system where both the entropy and energy are com- 
posed using non-additive rules 



3. Results 

When fitting data on fragmentation, we used the 
parametrisation of the TS distribution: 



1 do- 
er dx 



Ax ' 1 1 1 + 
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77(V5/2) 
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whith the parameters q — 1 + I /(a + D + 1) and 
T = ( yfs/2)J3/(D (a + D + 1) ). Fittings of Eq. @ to 



measured fragmentation functions (Fig.[T]l and ratios of 
data and fits (Fig. |2]i show that the one-dimensional TS 
reproduces data for about 0.01 ~ ^ ~ 0.2. The increase 
of q (Fig. [3} and decrease of T (Fig. |4j parameters with 
yfs is similar to what was found in ifTTl . and differs from 
pp results [16], where T was found to be independent 
of -\fs (though the analysis was done only for -\fs > 200 
GeV). 



Rewriting the microcanonical TS (Eq. ( 1 1 1) in terms 
of the parameters, q and T, we arrive at 



1 dcr 
o~ dx 



Ax°- ] (1 -x) 



i/W-i) ' 



(13) 



Fittings of this microcanonical generalisation of the TS 
in one-dimension to fragmentation functions (Fig. [SI as 
well as ratios of data and fits (Fig.|6]l, show that Eqf(13]l 
gives a good fit in a wide range, 0.01 J$ x 1. The in- 
crease of q with disappears (Fig.|7|, and the decrease 
of T with tJs reduces significantly (Fig. [8}. 

Both the canonical and the microcanonical TS in one 
dimension fail to describe data for very low x (x ;$ 
0.01). Introducing a scaling in the dimension parame- 
ter, D ~ 1 / yfs solves this problem, as can be seen in 



Figs. 17 and [18] showing fits. The interpretation of this 



is not straightforward, though. It should not mean that 
hadron production in a jet at -\fs = 200 GeV is a one- 
dimensional process, while at -\[s — 14 GeV it is a 
three-dimensional one. Let us consider that transverse 



momenta (pj) of hadrons in a jet (transverse with re- 
spect to the jet axis) are smaller than pro - 2 GeV/c. 
Then the ratio of the maximal transverse hadron mo- 
mentum and the maximal longitudinal hadron momen- 
tum (piQ = -\fs/2 considering 2-jet events) grows from 
Pto/Plo ~ 0.02 up to pro/pLO ~ 0.3 as yfs decreases 
from V* = 200 GeV down to yfs = 14 GeV. Hadrons 
with momenta p < pro (which means x < proli. V*/2)) 
may be produced isotropically in the hemi-sphere of the 
jet, while hadrons with momenta p » p T o are pre- 
destined to go approximately parallelly to the jet axis. 
So, the effective dimension of the phase space depends 
on p/ptq, rather than on x or on -\[s. Thus this causes 
scaling violation of the fragmentation functions. In 
non-perturbative models lt30Tl — If32]| a similar, (m T /Ej et ) a 
(with a ~ 1 — 2), energy dependence of spectra (mj 
being the transverse energy of the hadron, of order 0.5 
GeV) was reported as well. 

In Section |2j it has been shown that a Gamma- 
distribution of the shifted multiplicity N - No can result 
in a TS or microcanonical TS shaped spectrum. Fits of 
Eq. Q (s hown in Figs. 9p2 1 and of Eq. ( 10 1 (shown in 
Figs. 13p6] l to data on multiplicity distributions mea- 
sured at various collision energies, show that both as- 
sumptions, No = and jVo = 1 +2/D (with D — 1) give 
acceptable agreement with data. However, the predicted 
power of the spectrum from Eq. (Ill, scattering around 



12 and 10 (Figs. 1 1 and 15 I, is much higher than what is 



seen from fits to data on fragmentation functions. The 
fitted power, l/(q - 1) in Eq. ( 12 1 and in Eq. ( p"3j ), scat- 
ters around 2 and 2-4 respectively (Fig.[3]and Fig.|7]i. 



4. Conclusions 

In this paper, we compared inclusive data on e + e~ — > 
ft* + X reactions with the models on statistical hadron 
production described in Section [2] In these models, we 
considered 2-jet events, neglected the masses of quarks 
and hadrons. Furthermore we assumed that hadrons 
in a single jet are produced according to a statistical 
distribution. We used two types of distributions: the 
Boltzmann-Gibbs, and the microcanonical distribution 
in one dimension, for the reasons discussed in Section[2] 

We found that event-by-event fluctuations of the 
multiplicity can be described well by the Gamma- 
distribution, and that such fluctuations can result in 
an average fragmentation function of the form of the 
Tsallis-like, or a possible microcanonical generalisation 
of the Tsallis distribution. In Section [3] we have shown 
that Tsallis distribution reproduces data on fragmenta- 
tion for 0.01 ^ x ^ 0.2, while the microcanonical 
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Figure 1: Fragmentation functions measured at various collision 
energies (data of graphs from top to bottom are published in refs. 
Q3E1E3EZ1EIDE1ESIESIESIEI]) and fitted l dimensional TS 
distributions (Eq. \12\ with D = 1). 



Figure 3: Fitted values of the q parameter in Eq 
measured fragmentation functions shown in Fig.fT 



121 with D = 1 to 



If 




Figure 2: Ratios of measured fragmentation functions and fitted 1 
dimensional TS distributions (Eq. \ \2\ with D = 1) at various collision 
energies (data of graphs are published in refs. |.35 ]-[41 ]). 



Figure 4: Fitted values of the T parameter in Eq. j 1 2| with D = 1 to 
measured fragmentation functions shown in Fig.fT] 
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Figure 5: Fragmentation functions measured at various collision 

energies (data of graphs from top to bottom are published in refs. Fi § ure 7: Fitted values of me 1 parameter in E q. (13) with D = 1 to 

Q3|lJ|33|32]|3l]|Il|I3|I3|I3III]) and fitted 1 dimensional MTS measured fragmentation functions shown in Fig.|3] 
distributions (Eq. {l3j with D = 1). 
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Figure 6: Ratios of measured fragmentation functions and fitted 1 
dimensional MTS distributions (Eq. (13) with D = 1) at various colli- 
sion energies (data of graphs are published in refs. |35|-|41 ]). 



Figure 8: Fitted values of the T parameter in Eq. (13) with D = 1 to 
measured fragmentation functions shown in Fig. [3] 
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Figure 9: Multiplicity distributions measured at various collision en- 
ergies (data of graphs from top to bottom are published in refs. I42I 
|42]|42j|42j|43]|44]|34][44]|44]|^|44j^ ) 
and fitted Gamma-distributions (Eq. 4 1 0| with Nq = 0). 



Figure 1 1 : Fitted values of the a parameter in Eq. (Tolwith Nq = 
to the measured multiplicity distributions shown in Fig. 191 
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Figure 10: Ratios of measured multiplicity distributions and fitted 
Gamma-distributions (Eq. \10\ with Nq = 0) measured at various col- 
lision energies (data of graphs are published in refs. |40| and |42|- 
ED). 



Figure 12: Fitted values of the /J parameter in Eq. j 1 0| w ith Nq = 
to the measured multiplicity distributions shown in Fig. 191 
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Figure 13: Multiplicity distributions measured at various collision 
energies (data of graphs from top to bottom are published in refs. [42 

|42j|42j|ia|43l|44]|44][44]|44]|^|44j|44^ 

and fitted shifted Gamma-distributions (Eq. \IQ) with Nq = 3). 



Figure 15: Fitted values of the a parameter in Eq. (Tolwith No 
to the measured multiplicity distributions shown in Fig.|13| 




Figure 14: Ratios of measured multiplicity distributions and fitted 
shifted Gamma-distributions (Eq. \10\ with JV U = 3) measured at var- 
ious collision energies (data of graphs are published in refs. [40 1 and 

El-GOD). 



Figure 16: Fitted values of the /J parameter in Eq. j 1 0| wi th Nq = 3 
to the measured multiplicity distributions shown in Fig.|13| 
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Tsallis for 0.01 



< 



< 



1. Below 



0.01, the one- 



dimensional models tested, deviate from measurements, 
especially for -\fs < Mz- 

Introduction of an effective phase-space dimension, 
scaling like D ~ 1 / -\fs seemingly solves the prob- 
lem. However, the underlying physics, as argued in Sec- 
tion [3] supports the emergence of a scale-breaking term 
in the fragmentation function of the form, (mjjEj et ) a , as 
reported in Refs. [30|-[32|, rather than a scaling of the 
dimensionality of the phase space. Furthermore, when 
considering low-x hadron production, yields from 3-jet 
events are not negligble, especially for yfs < Mz- 

Anyway, the decisive argument for or against the 
models introduced above would be the measurement 
of momentum distribution of hadrons in single jets for 
fixed multiplicity. 
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